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We discuss a generalization of the surface reconstruction arguments for the struc- 
ture of intermediate size Si clusters, which leads to model geometries for the sizes 
33, 39 (two isomers), 45 (two isomers), 49 (two isomers), 57 and 61 (two isomers). 
The common feature in all these models is a structure that closely resembles the 
most stable reconstruction of Si surfaces, surrounding a core of bulk-like tetrahe- 
drally bonded atoms. We investigate the energetics and the electronic structure of 
these models through first-principles density functional theory calculations. These 
models may be useful in understanding experimental results on the reactivity of Si 
clusters and their shape as inferred from mobility measurements. 



I. INTRODUCTION 



Silicon is a material that attracts consider- 
able interest due to its technological importance. 
It has also come to be regarded as the repre- 
sentative covalently bonded solid. In the last 
decade, a new form of Si has generated much 
excitement because it holds promise for novel 
technological applications as well as for provid- 
ing unique insight, into the nature of covalently 
bonded materialstl. This new form is clusters of 
Si, which consist of a few tens to a few hundreds 
of atoms and have properties different from the 
bulk. The properties of Si clusters depend on 
their size, and a detailed study of this depen- 
dence could potentially lead to a better under- 
standing of how chemical bonding evolves from 
that characteristic of small molecules to that 
characteristic of the bulk. 

Despite much theoretical and experimental 
work on Si clusters during the last decade, little 
is known about their structure. It is self evident 
that for very small cluster sizes, all of the atoms 
in a cluster will be exterior atoms, in the sense 
that there is no shell formed by a subset of atoms 
that completely surrounds any one atom. These 
very small Si clusters should therefore be viewed 
as molecules, i.e. entities with unique structural, 
physical and chemical characteristics that do not 
resemble other forms of Si. TJiis picture was 
originally proposed by PhillipsEI, who analyz 
early experiments on fragmentation patternsu 
Equally evident is the fact that beyond a cer- 
tain size there will be interior atoms completely 
surrounded by a shell of other exterior atoms. 
(In the following, we use the terms "interior" 
and "exterior" as defined above to characterize 



cluster atoms, and reserve the terms "bulk" and 
"surface" to characterize atoms of the solid.) It 
is also natural to expect that for a large enough 
size, the cluster will resemble a crystalline frag- 
ment. 

The size at which the transition from 
molecules to bulk-like fragments takes place has 
not yet been determined precisely, although se 
eral theoretical predictions have been madcfl' 
It has also been proposed that the experimen- 
tally observed sharp transition in the shape of 
Si clustersQ, which occurs at around size 27 may 
be related to the onset of structures with in- 
terior ("bulk" -like) and exterior ("surface" -like) 
atoms in the clustem. Such a transition in 
shape had actually been predicted by theoret- 
ical simulationdj. 

The structure of the small Si clusters (up 
to 10 atoms) has been deterpuued by ex- 
tensive theoretical calculationgliH~t3 and by 
comparison cf£, calculations to spectroscopic 
measuremental3~EII. Much less is known for the 
larger clusters, although a number of interest- 
ing models have been proposed. One po ^ibiL - 
ity is models analogous to fuUerene cagestSEj, 
but given the tendency of Si to form four-fold 
coordinated structures this seems unlikely. A 
different approach is to consider models. [With 
non-spherical shapes for medium sizesEl£3'Eil or 
withjiiterior and exterior atoms for larger sizes, 
butE3"E3. The assumption that there exist in- 
terior and exterior atoms does not necessarily 
imply that their respective environment will re- 
semble that of bulk or surface atoms. Never- 
theless, it is an attractive proposition to con- 
sider what types of structures can be produced 
if geometries resembling the bulk structure (for 
interior cluster atoms) and the surface recon- 
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structions (for exterior cluster atoms) were to 
dominate the cluster geometry. Only a few 
cluster sizes are compatible with the require- 
ment that all interior cluster atoms have envi- 
ronments that closely resemble the bulk struc- 
ture, and all exterior cluster atoms have envi- 
ronments that closely resemble surface recon- 
structions. We call the resulting models Surface 
Reconstruction Induced Geometries (SRIG's for 
short). This idea has been invoked as an ex- 
planation |Of the existence of magic numbers in 
Si clustera23 as revealed-bjj-their reactivity with 
various chemical agent£jE3. 

In this paper, we present in Section || a thor- 
ough investigation of SRIG models for the struc- 
ture of selected clusters spanni ng the sizes from 



33 to 61 atoms. In Section III we use first- 
principles density functional theory calculations 
to study the energetics and electronic proper- 
ties of these models. Finally, in Section |^ we 
comment on the relevance of these models to 
experimental observations of the reactivity and 
the mobihty of Si clusters. 



II. SURFACE RECONSTRUCTION 
INDUCED GEOMETRIES 

Based on the assumption that the surface re- 
construction determines the cluster geometry, 
we have constructed models for Si clusters in 
which exterior atoms resemble closely the coor- 
dination and bonding of atoms in the most sta- 
ble surface reconstructions of Si, including the 
(111) 7x7 diper-adatom-stacking fault (DAS) 
reconstructiorH, the (100) 2x1 dimer recon- 
struction, andj-^e (111) 2x1 7r-bonded chain 
reconstructiorH. All interior atoms have the 
coordination and bonding arrangement of bulk 
Si atoms, i.e. four-fold tetrahedral coordination. 
The special sizes for which we were able to con- 
struct such models include the sizes 33, 39, 45, 
49, 57 and 61. This does not exclude the pos- 
sibility of other sizes that have similar bonding 
arrangements in this size range. In order to con- 
struct these models we imposed two additional 
constraints: (i) that the clusters have an atom at 
their geometric center, and (ii) that the clusters 
have tetrahedral overall symmetry. These con- 
straints were motivated by simplicity considera- 
tions and not from any physical requirements. If 
the constraints are relaxed, it might be possible 
to construct additional structures with similar 



characteristics but lower symmetry. 

In describing the cluster geometries, we adopt 
the following conventions: We will display only 
the exterior atoms of the cluster, in a two- 
dimensional picture which is obtained by un- 
folding three of the four sides of the tetrahe- 
dron symmetrically around the fourth side. In 
so doing, one obtains an equilateral triangle, at 
the center of which there is a smaller, inverted 
equilateral triangle corresponding to one side of 
the tetrahedron around which the other three 
sides have been unfolded. We denote the sides of 
the tetrahedron by dotted lines in all the struc- 
tural figures. In the models we discuss there 
exist one to five interior atoms, which are not 
shown in the structural figures. The position of 
these interior atoms is discussed in detail along 
with every structural model. What is actually 
shown in the structural figures are the projec- 
tions of atoms on the faces of the tetrahedron. 
For simplicity we refer to these projections as 
the atomic positions. This manner of displaying 
the cluster geometries has the added advantage 
that the reader can easily construct schematic 
three-dimensional models of the cluster exterior 
by cutting out the two-dimensional figures, fold- 
ing them at the three edges of the central in- 
verted triangle, and gluing the outer edges to- 
gether. 

We have also adopted a shading scheme which 
is common to all the cluster sizes: First, we indi- 
cate by open circles (we call those white atoms) 
the exterior atoms that are at the apexes of the 
tetrahedron; there are four such atoms in each 
model. In some of the models (45B, 61A, 61B) 
there exists an additional set of four atoms at the 
centers of the faces of the tetrahedron, which we 
also denote by open circles. Next, we indicate 
by hatched circles the atoms that are bonded to 
the first set of white atoms; there are twelve such 
atoms in each model. Next, there are atoms that 
are bonded to the twelve hatched atoms, and 
are indicated by grey circles. Finally, there are 
atoms shown in black, which are bonded to the 
grey atoms or to the white face-center atoms, 
but not to the hatched or white apex atoms. 
Sets of atoms with the same shading are equiv- 
alent by symmetry. With these conventions, we 
proceed to describe the individual cluster mod- 
els in detail. 

We begin with the cluster of size 33, shown 
in Fig. 1_ This cluster had been proposed 
previouslyE3, and consists of 28 exterior atoms 
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and 5 interior atoms. Of the five interior atoms, 
one is at the geometric center of the tetrahe- 
dron and the other four are bonded to the cen- 
tral atom in directions pointing away from the 
center toward the vertices of the tetrahedron. 
An interior atoms are four-fold coordinated as 
atoms in bulk Si are. The neighbors of the four 
interior atoms other than the central one are the 
twelve hatched exterior atoms. The hatched ex- 
terior atoms are also four-fold coordinated, but 
their bonds are at angles severely distorted rel- 
ative to the tetrahedral angle. The neighbors 
of the hatched exterior atoms are the four in- 
terior atoms (other than the central one) and 
the twelve grey exterior atoms which are three- 
fold coordinated: they have two bonds to the 
hatched exterior atoms and one bond to another 
grey exterior atom each. Finally, there are four 
exterior white atoms at the apexes of the tetra- 
hedron; these atoms are three-fold coordinated, 
and are bonded to three hatched exterior atoms. 

The close similarity of this structure with sur- 
face reconstructions of Si involves two elements: 
(a) The four white apex atoms are in local ge- 
ometries which are very similar to the adatoms 
on the Si(lll) 7x7 reconstruction: the adatoms 
are directly above second layer bulk Si atoms (in 
the present model this role is played by the four 
interior atoms other than the central one), (b) 
The twelve grey exterior atoms are bonded in 
pairs like the dimers on the Si(lOO) 2 x 1 re- 
construction; these atoms also play the role of 
the dimers that characterize the Si(lll) 7x7 
reconstruction and help stabilize the adatom re- 
construction. Thus, in this model there is a co- 
operative effect that combines features from the 
two most stable surface reconstructions of Si to 
produce a highly symmetric structure with 5 in- 
terior and 28 exterior atoms. 

It is noteworthy that the bonding arrange- 
ment of exterior atoms in this structure consists 
of pentagonal rings (with corners at one white, 
two hatched and two grey atoms) and distorted 
hexagonal rings (they appear to be almost trian- 
gles in Fig. 1, with corners at three hatched and 
three grey atoms). In this sense, the shell of ex- 
terior atoms resembles the structure of the C28 
fuUerene, which has been studied-, extensively as 
the smallest stable fullereneE3~E3. The resem- 
blance, however, is superficial, since it is limited 
to this geometric aspect, while the C28 fuUerene 
is stabilized by the 7r-bonding of atoms at the six 
hexagonal rings and the 8133 structure is stabi- 



lized by the presence of the adatom and dimer 
surface features, as well as the presence of the 
five interior atoms. Additional important dif- 
ferences concern the electronic properties of the 
two models, wiih the C28 cluster being an open 
shell structureEj, while the 8133 model is a closed 
shell structure (see Section |^) . 

We next consider the two models for the Si 
cluster with 39 atoms. The first model (referred 
to as 39A, shown in Fig. 2) is similar to the 
33 atom model discussed above, with six addi- 
tional atoms that form bonds to the grey atoms 
(shown as black, at the midpoints of tetrahe- 
dron edges). In this model, the white, hatched 
and grey atoms are coordinated as in the 33 
atom model. The six extra atoms are placed be- 
tween pairs of grey atoms, and form new bonds 
to those atoms which become four-fold coordi- 
nated. The six extra atoms are two-fold coor- 
dinated each. This type of atomic arrangement 
is not usually observed on Si surfaces. How- 
ever, it is a reasonable arrangement for Si atoms 
which cannot find close neighbors to bond to. 
The preferred bonding for these two-fold coor- 
dinated atoms will be through p orbitals to their 
two neighbors, while they retain two of their four 
valence electrons in a low-energy non-bonding s 
orbital. 

The second model for the 39 atom cluster (re- 
ferred to as 39B, shown in Fig. 3), consists of 
a different geometric pattern. The four white 
apex atoms are in similar arrangements as in 
the 39A model. The hatched atoms have three 
bonds, one to interior atoms, one to white apex 
atoms and one to black atoms at the centers 
of the tetrahedron edges. The grey atoms no 
longer form dimers; instead, they form trimers 
centered at the faces of the tetrahedron. This is 
a somewhat unusual arrangement for Si atoms, 
and is not encountered in native surface recon- 
structions of Si, since the presence of trimer 
units induces significant strain on the surface. 
However, a trimer reconstruction of this type 
is common for group- \L-atoms on the surfaces 
of III-V semiconductorsLj, or when the group-V 
elements are u^d as passivating layers on the 
Si(lll) surfaceEII. Moreover, group-IV elements 
like Sn and Pby-ialso form trimer units on the 
Si(lll) surfaceEa. We suggest that under the 
proper conditions (represented here by the size 
of the cluster) similar trimer units of Si atoms 
may be stable on the cluster surface. Finally, 
there are six atoms at the centers of the edges 
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of the tetrahedron, which in this model arc four- 
fold coordinated with two bonds to the hatched 
atoms and two more bonds to grey atoms in the 
trimers. This arrangement results in three-fold 
and seven-fold rings on the exterior shell of the 
cluster and is markedly different from the previ- 
ous geometry for the size 39. 

The first model we considered for the 45 atom 
cluster (referred to as 45A, shown in Fig. 4) 
is a simple variation of the 39B model: the 
trimers at the centers of the tetrahedron faces 
are replaced by hexagons and the black atoms 
at the centers of the tetrahedron edges are elim- 
inated. This should help reduce the strain in- 
duced by the trimers, and it should enhance the 
stability of the apex atoms: the coordination 
of their neighboring hatched atoms has been re- 
stored to four (one apex atom, one interior atom 
and two grey atoms), while the exterior pattern 
again contains three pentagonal rings around 
each white apex atom and dimers of grey atoms. 
In addition to the pentagonal rings (with corners 
at one white apex, two hatched and two grey 
atoms), the exterior of the 45A model contains 
two types of hexagons: one composed of six grey 
atoms and one composed of four grey and two 
hatched atoms. In this sense, the exterior of this 
model bears a supesficial resemblance to the C40 
fuUerene structura2il, but just like in the case of 
Si33 and its companion fuUerene C28, the stabil- 
ity and electronic properties of the two types of 
clusters are quite different. 

The next model for a 45 atom Si cluster (45B) 
is illustrated in Fig. 5, and corresponds to a dif- 
ferent type of exterior bonding arrangementE3. 
In this case there is only one interior atom, 
bonded to the four white atoms that occupy 
the centers of the tetrahedron faces. The white 
atoms at the tetrahedron apexes are no longer 
at positions that resemble the adatom struc- 
ture of the Si(lll) 7x7 reconstruction. In- 
stead, they form intersection points where zig- 
zag chains consisting of hatched and grey atoms 
meet. The presence of these chains should pro- 
duce a low energy exterior shell because they re- 
semble closely the chains of atoms in the other 
stable reconstruction of the Si(lll) surface, the 
TT-bonded chain model proposed by PandeyE3. 
In this model, the dangling bonds on the neigh- 
boring chain atoms form 7r-bonded combina- 
tions. The geometry of the exterior shell con- 
sists of pentagonal rings around the four white 
apex atoms and hexagonal rings around the four 



white face-center atoms. Since there is only 
one interior atom bonded to the four face-center 
white atoms, the hatched atoms are only three- 
fold coordinated, with all their neighbors be- 
longing to the exterior shell. The grey and 
black atoms are three-fold coordinated, as are 
the apex white atoms. The only atoms that 
are four-fold coordinated, other than the central 
atom, are the four white atoms at the centers of 
the tetrahedron faces. This model has chirality 
(i.e. there exist equivalent left and right handed 
versions), which is evident from the fact that the 
zig-zag chains do no lie on any high symmetry 
direction of the tetrahedron, and thus break the 
left-right symmetry of the previous models. 

The same exterior shell can be used to con- 
struct two different models for the 49 atom clus- 
ter. We refer to these models as 49A and 49B. 
The first (49A) constists of the exterior shell 
shown in Fig. 5 plus five interior atoms, one at 
the geometric center of the tetrahedron and four 
more bonded to the central atom, and point- 
ing away from it toward the four white apex 
atoms. The second (49B) consists of the same 
exterior shell and five interior atoms again, one 
at the geometric center of the tetrahedron and 
four more bonded to the central one and point- 
ing away from it, toward the four white atoms 
at the centers of the tetrahedron faces. In the 
case of model 49A, the four white apex atoms 
have a coordination similar to that of adatoms 
on the Si(lll) 7 X 7 DAS reconstruction, as de- 
scribed in the case of the 33 atom model. In 
the case of model 49B, the role of adatoms is 
assumed by the four white atoms at the cen- 
ters of the tetrahedron faces. In both cases, the 
atoms that are bonded to the adatoms become 
four-fold coordinated. Specifically, in model 
49A the hatched atoms become four-fold coordi- 
nated, while in model 49B the black atoms be- 
come four-fold coordinated. Since the exterior 
geometry in models 49A and 49B is the same as 
in the model 45B, the pattern of pentagons and 
hexagons is the same. The important difference 
between models 49A and 49B is that in the first 
the adatoms (white apex atoms) are surrounded 
by three pentagonal rings, while in the second 
the adatoms (white face-center atoms) are sur- 
rounded by three hexagonal rings. 

The model for the 57 atom cluster is an ex- 
tension of the 45A model, in which the six-fold 
rings of grey atoms on the tetrahedron faces are 
replaced by nine-fold rings as shown in Fig. 6. 
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This is achieved by the addition of three ex- 
tra atoms (shown in black) on each tetrahedron 
face. The rest of the cluster structure is the 
same as the 45A model, with five interior atoms 
(one at the geometric center of the tetrahedron 
and four more bonded to it, pointing toward the 
four white apex atoms) and the four white apex 
and twelve hatched atoms in similar positions 
as before. The extra atoms that turn the face- 
centered hexagonal rings of the 45A model into 
nine-fold rings, are themselves bonded in dimer- 
like pairs across the the centers of the tetrahe- 
dron edges. In this way the exterior shell is com- 
posed of nine-fold rings centered at the tetrahe- 
dron faces and two types of pentagonal rings, 
one with corners at a white apex, two hatched 
and two grey atoms, and one with corners at a 
hatched, two grey and two black atoms. 

The two models we have constructed for the 
61 atom cluster involve geometries that have 
four white face-center atoms, in two different 
configurations. The first configuration (labeled 
61A, shown in Fig. 7) derives from the 57 
atom model, with the additional white atoms 
bonded to the three black atoms on each face, 
and with the dimer bonds between black atoms 
broken. The remaining atoms are in exactly the 
same configuration as in the 57 atom model. 
This change modifies the exterior shell, which 
now consists of eight-fold rings (centered at the 
edges of the tetrahedron and composed of two 
hatched, four grey and two black atoms) and 
two types of pentagonal rings, one composed of 
a white apex, two hatched and two grey atoms, 
and the other composed of a white face-center, 
two black and two grey atoms. In this model, as 
far as the exterior shell is concerned, the white 
apex atoms and the white face-center atoms are 
equivalent. Therefore the four interior atoms 
(other than the central one) can be considered 
to be pointing toward either the white apex or 
the white face-center atoms. 

Finally, the second model for a 61 atom clus- 
ter (labeled 61B, shown in Fig. 8) contains 
four white apex atoms at adatom positions, sur- 
rounded by hatched atoms that form pentag- 
onal rings with adjacent grey atoms, and four 
white face-center atoms that form hexagonal 
rings with a hatched, two grey and two black 
atoms. In this case the four white apex atoms 
and the four white face-center atoms are not 
equivalent. The four interior atoms (other than 
the central one) are bonded to the four white 



apex atoms, since those have the benefit of being 
surrounded by pentagonal rings which should 
stabilize their adatom-like features. In the 61B 
model, the exterior shell consists of pentagonal 
and hexagonal rings, which gives it a superficial 
resemblance to the corresponding C56 fuUerene 
structure. 



III. STRUCTURAL, ENERGETIC AND 
ELECTRONIC FEATURES OF SRIG 
MODELS 

Having described the geometric features of the 
SRIG models, we discuss next the results of en- 
ergy optimization and electronic structure cal- 
culations for these models. The calculations are 
based on density functional theory in the local 
density approximatiorO (DFT-LDA. Although 
these calculations do not necessarily provide the 
most accurate comparisons of cluster energies, 
they are reasonably reliable in determining opti- 
mal geometries by minimizing the magnitude of 
the calculated Hellmann-Feynman forces. The 
present calculations employ a plane wave ba- 
sis with a cutoff energy of 8 Ry, and a cubic 
supercell with lattice constant equal to 15.875 
A(which gives a basis of 10400 plane waves), 
for the clusters of sizes 33 to 49, and equal to 
17.463 A, for the clusters of sizes 57 and 61 
(which gives a basis of 13600 plane waves). With 
these supercells the periodic images of the clus- 
ters are reasonably well separated. A single k- 
point (the center of the supercell Brillouin Zone) 
was used in all the cluster calpjlations. The 
Car-Parrinello iterative schemeCj was used to 
solve the Kohn-Sham equations self consistently, 
and non-local norm conserving pseudopotpitials 
were employed to model the atomic coreaij. 

In Table I we provide a list of the atomic po- 
sitions for the SRIG models. In this Table only 
the positions of inequivalent atoms are given, 
with the remaining positions obtained by ap- 
plying symmetry operations of the tetrahedral 
group. The position of the central atom in each 
cluster is taken to be the origin of the coordinate 
system and is omitted from the list of positions. 

We discuss first certain structural features of 
the optimized cluster geometries. These con- 
siderations were motivated by some exciting re- 
cent experiments that are able to resolve as- 
pects of the cluster geometry through gas phase 
mobility measurementsc3. In the simplest pic- 
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ture, the mobility measurements reveal the pro- 
jected cross section of the cluster. For spheri- 
cally shaped clusters, this cross section is given 
by = 7r6^j„, where bmin is the distance of clos- 
est approach between gas phase molecules and 
the cluster. By considering the atoms on the 
surface cluster as points, we can obtain 6„im as 
the distance of the atom farthest from the clus- 
ter center, which gives the values of f2 tabulated 
in Table II. This quantity is a coarse measure 
of the cluster geometry and gives approximately 
the overall cluster size in cross section. From the 
comparison of O values in Table II, it is evident 
that the cluster cross section does not follow ex- 
actly the number of atoms in the cluster. It is 
also interesting that the clusters with the high- 
est cohesive energy (see following paragraph and 
Table II) also have the smallest (especially the 
models 33, 45 A, 49A and 49B), i.e. the lowest 
energy clusters are also the most compact ones. 

As, Shvartsburg and Jarrold have pointed 
ouio, the projection of the cluster cross sec- 
tion is appropriate for clusters which are locally 
convex. When the surface of the cluster con- 
tains concave regions, the scattering cross sec- 
tion is increased due to multiple collisions. In 
order to obtain a measure of this effect, we have 
calculated the ratio Rc of surface triangles that 
are concave to the total number of triangles on 
the cluster surface. Since this is a measure of 
local concavity, we have limited the definition 
of surface triangles to those which have sides 
smaller than 6.35 A. Table II gives the values 
of Rc for the various models. Two of the clus- 
ters, 33 and 45A, are everywhere locally con- 
vex {Rc — 0). For the rest of the clusters, Rc 
ranges between 1.9% to 10.8% which is relatively 
small, but not insignificani, since the accuracy 
of mobility measurementsCJ is about 2%. The 
clusters with the smallest non-zero Rc are 49A 
and 49B, which together with the clusters 33 
and 45A form the group of clusters which are 
most compact (lowest fi) and have the highest 
cohesive energy. These comparisons should be 
helpful in identifying which of the SRIG models 
may be relevant to the structure of Si clusters 
as determined by mobility experiments. 

For the energetic comparisons, since the type 
of calculations employed here do not provide an 
accurate estimate of the cohesive energy of bulk 
Si, we have opted to quote the cohesive energy 
per atom of each cluster by comparing the clus- 
ter energy to an equivalent number of bulk Si 



atoms, and using the experimental value for the 
cohesive energy of the bulk (4.68 eV). In this 
way, the quoted energies per atom for the clus- 
ters should be closer to realistic cohesive ener- 
gies. For the bulk calculation, we use the same 
plane wave cutoff, the primitive 2-atom unit cell 
of the diamond lattice, and a set of k-points that 
produces a density in reciprocal space equivalent 
to the center of the supercell cube in the cluster 
calculations. 

Table II contains a comparison of the cohesive 
energy per atom and the corresponding energy 
gap between the highest occupied molecular or- 
bital (HOMO) and lowest unoccupied molecu- 
lar orbital (LUMO), for the different clusters we 
studied. The model with the highest cohesive 
energy per atom (49A) has almost 82% of the 
bulk cohesive energy, while the model with the 
lowest cohesive energy (39B) has 77% of the bulk 
cohesive energy. Fig. 9 shows the density of 
states (DOS) for each model. The symmetry of 
the clusters dictates that the electronic states 
are singly, doubly or triply degenerate, accord- 
ing to the dimensions of the irreducible repre- 
sentations of the tetrahedral group. In all cases 
the total band width (between 12 and 13 eV) is 
comparable to that of the bulk. 

There are some interesting insights revealed 
by these comparisons of cohesive energies and 
electronic structure. The cluster with the small- 
est cohesive energy (least stable energetically) is 
the cluster with 39 atoms and trimer units at the 
tetrahedron surfaces (model 39B, Fig. 3). The 
relatively high energy of this cluster can be ra- 
tionalized as being due to the strain induced by 
the trimer units. Interestingly, this cluster also 
has the largest HOMO-LUMO gap. Thus, at 
least for the structures considered here, there is 
no direct correlation between the energetic sta- 
bility and the HOMO-LUMO gap. The exis- 
tence of a large HOMO-LUMO gap may be in- 
dicative of low chemical reactivity, since in this 
case the system has a closed electronic shell. In 
this sense, if the 39B SRIG model were to be 
formed, it might be expected to have low chem- 
ical reactivity. However, a quantitative study 
of chemical reactivity should involve a more de- 
tailed examination of cluster electronic states, 
as well as case studies of chemical reactions be- 
tween |4ie cluster and representative molecular 
ag entsB The presence of a different model of 
the same size (39A) which has lower energy and 
a comparable HOMO-LUMO gap suggests that 
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the hkehhood of the geometry 39B being real- 
ized is small. 

Another interesting point is that the low- 
est energy clusters are the two models with 49 
atoms, closely followed by the 45A model and 
the 33 model. While the 45A and 33 models 
have significant HOMO-LUMO gaps, the two 
49-atom models have very small (49A) or non- 
existent (49B) gaps. Thus, based on the sim- 
ple arguments mentioned above on the relation 
of the HOMO-LUMO gap to chemical reactiv- 
ity, one might expect that the models 33 and 
45A will exhibit low chemical reactivity, while 
models 49A and 49B will have higher chemical 
reactivity, despite their lower energy per atom. 
This is consistent with the experimental find- 
ings of Jarrold and Honeao, who observed that 
chemical reactivity and thermodynamic stabil- 
ity are not related. Finally, the larger clusters 
of sizes 57 and 61 again have very small (57) or 
non-existent (61A, 61B) gaps, suggesting high 
chemical reactivity despite their relative ener- 
getic stability (especially for model 61 A). 



IV. DISCUSSION AND CONCLUSIONS 

We have presented a detailed discussion of the 
geometric features, the relative energies and the 
electronic structure of models for Si clusters of 
sizes 33, 39, 45, 57 and 61. The basic char- 
acteristic of these models is their resemblance 
to bulk Si, including geometries of the exte- 
rior atoms that resemble surface reconstruction 
features and geometries of the interior atoms 
that resemble the four-fold tetrahedral coordi- 
nation of the bulk. Although this is an appeal- 
ing feature, it does not guarantee that these 
models are the energetically preferred geome- 
tries. In fact, it seems that molecular dynam- 
ics (MD) or simulated jannealing (SA) simula- 
tions based on empiricalo, .sftpi-empiricaO and 
first-principles calculationsOEJ tend to give ge- 
ometries different from the ones proposed here. 
Therefore, it is important to address to what ex- 
tend the models considered here are relevant to 
the structure of real Si clusters. 

We suggest that these models may indeed be 
relevant to the structure of real Si clusters for 
the following reasons: 

(1) It appears from experiment that there 
is something special about the chemical re- 
activity of the (so called "magic") sizes 33, 



39 and 45, while other clusters have approxi- 
mately constant and much higher (by several or- 
ders of rnaenitude) reactivity from these magic 
numbersc3 EErEJ. It is appealing that for all the 
magic number sizes we were able to construct 
SRIG models, and that these particular mod- 
els have the largest HOMO-LUMO gaps of all 
the structures considered. Moreover, the SRIG 
models of sizes 49, 57 and 61, which are energeti- 
cally equally stable stable as the magic numbers, 
have very small or non-existent HOMO-LUMO 
gaps. In this sense, the pattern of SRIG mod- 
els is compatible with the pattern of experimen- 
tal measurements of reactivity, assuming that 
the HOMO-LUMO gap can be used as a coarse 
measure of chemical reactivity (with the caveats 
mentioned in Section ^ ) . 

(2) If the clusters of various sizes possess no 
special geometric features, as the simulatioas 
based on various methodologies suggestc^E^, 
then it is difficult to explain the dramatic 
changes in reactivity and the existence of the 
magic numbers. Actually, MD simulations 
based on DFT-LDAc3, reported that by aug- 
menting by one atom the core of their opti- 
mal 45-atom geometry consisting of a 38-atom 
outer shell and a 7-atom core, a cluster of 46 Si 
atoms is obtained which has the same cohesive 
energy, the same outer shell and the same num- 
ber of "danging bonds" as the original 45-atom 
model. In other words, adding one atom to the 
45-atom model gives a cluster of the same sta- 
bility (cohesive energy) , surface structure (outer 
shell) and electronic states (number of dangling 
bonds). Consequently, the reactivity of the 45 
and 46 atom clusters obtained by these simu- 
lations should be essentiallyi-idetttical, in direct 
contradiction to experimentL3~Ot^l. The SRIG 
models presented here give a natural explana- 
tion to this problem: since the SRIG models 
(e.g. the 45-atom model with the lowest energy) 
correspond to a perfectly reconstructed outer 
shell, the addition or subtraction of one atom 
will drastically change the structure because it 
cannot be accommodated as part of the outer 
shell reconstruction or of the interior, and thus 
the reactivity of the cluster will be significantly 
increased. 

(3) The methodologies employed to search for 
possible cluster structures in MD or SA simula- 
tions have limitations, especially when it comes 
to comparing geometries that are very differ- 
ent in bonding arrangements. This applies to 
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the empirical and semi-empirical methodologies, 
but also to a certain extend to the DFT-LDA 
methodologies. It has been argued, for exam- 
ple, that electron correlation effects are cru- 
cial in jSlabilizing the structure of certain Si 
clustersLJ'La, while such effects are not esnlicitly 
included in the DFT-LDA simulationsHEl. In 
the present work, we have used the DFT-LDA 
approach only for optimization of the structure 
of models within the proposed very strict ge- 
ometric constraints, that is, as if the overall 
structure of the cluster were known. It is well 
established that this methodology gives reason- 
able agreement with experiment for the detailed 
structural features (i.e. bond lengths and bond 
angles) of molecules, when the overall geometry 
is known. The electronic structure correspond- 
ing to a known geometry is also reasonably well 
obtained. There are however examples where 
this methodology has proven inadequate in com- 
paring the energetics of very different structures 
for a cluster of a given size and composition, for 
instance the cases,|J5there strong correlation ef- 
fects are presentE^Eil. In this sense, it may be 
overly optimistic to rely on DFT-LDA simula- 
tions in order to determine the structure of Si 
clusters in the size range considered here, with- 
out imposing any geometric constraints. If this 
is true, the results of empirical or semi-empirical 
simulations could be even less reliable, as far as 
unconstrained structure optimization of specific 
cluster sizes is concerned. 

(4) Even if a methodology that is very ac- 
curate were available, the task of determining 
the optimal structure of clusters in this size 
range is daunting. The review of Jones and 
GunnarsonE3 presents persuasive arguments on 
the intractability of finding the global energy 
minimum for a structure with several tens of 
atoms through an unguided search over configu- 
rational space. In fact, it is very likely that the 
structures with special properties (such as the 
magic numbers) correspond to deep and narrow 
wells in the multi-dimensional configurational 
space, which are difficult to locate by unguided 
searches. What we have provided here are phys- 
ically motivated models that could potentially 
correspond to such deep and narrow wells in the 
energy landscape. The proper approach would 
then be to use such models as starting points for 
searches of the configurational space. Even in 
this case, extreme caution should be used in per- 
forming simulations, since it is relatively easy to 



bring the structure outside the well if the bottom 
has not been reached by thorough relaxation. As 
an example, in the present work we have found 
that the cluster geometry can be stuck in local 
energy minima and is prevented from reaching 
the bottom of the energy well, simply due to 
the initial occupation of the electronic states. 
For this reason, the number of electronic states 
used in a DFT-LDA simulation must exceed the 
number of states required to accommodate the 
electrons, and the filling of the states near the 
Fermi level must be varied in order to achieve 
good relaxation and to avoid getting trapped in 
local energy minima. Once outside the mini- 
mum energy well, the simulation may end up 
exploring irrelevant and uninteresting structures 
which correspond to broad shallow basins in the 
energy landscape. 

In conclusion, we have discussed a set of mod- 
els for Si clusters in the range 33 - 61 atoms, 
which are physically motivated and exhibit in- 
teresting patterns in their energy and electronic 
structure. These models may be relevant to un- 
derstanding the exceptionally low chemical reac- 
tivity of certain magic number sizes. The pat- 
tern of HOMO-LUMO gaps exhibited by these 
models is consistent with the magic number 
clusters of low reactivity. Finally, the models 
are compact clusters with little concave surface 
area, and the lowest energy ones tend to have 
small cross sections. It will be interesting to 
consider how these models compare to the inna- 
vative experiments of Jarrold and coworkergl3, 
on the gas phase mobility of Si clusters. 
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TABLE I. Atomic positions in a.u. of representative atoms for the optimized structures of clusters of 
sizes 33 - 61. Atoms labeled 1 are in the interior of the cluster, the rest of the atoms are on the exterior 

and correspond to the atoms shown in Fig. 1 8. 





33 


39A 


39B 


45A 


45B 


49A 


49B 


57 


61A 


61B 


l{x) 


2.39471 


2.38887 


2.59713 


2.57556 


2.50590 


2.48468 


-2.84402 


2.47912 


2.64227 


2.52290 


i(y) 


2.39471 


2.38887 


2.59713 


2.57556 


2.50590 


2.48468 


-2.84402 


2.47912 


2.64227 


2.52290 




2.39471 


2.38887 


2.59713 


2.57556 


2.50590 


2.48468 


-2.84402 


2.47912 


2.64227 


2.52290 


2{x) 


5.21249 


5.16076 


2.47377 


2.99354 


-0.11844 


0.85354 


1.15471 


3.09986 


3.16035 


1.74224 


2(2/) 


5.21249 


5.16076 


2.47377 


2.99354 


3.90217 


4.39152 


4.53560 


3.09986 


3.16035 


5.59316 


2{z) 


0.25158 


0.23991 


-5.51192 


7.03966 


5.80660 


6.28003 


6.21782 


6.85742 


7.11649 


5.59316 


3(x) 


1.52358 


1.65023 


3.39552 


1.18201 


1.73935 


-2.03724 


-1.72237 


1.10146 


1.13596 


1.31967 


3(y) 


1.52358 


1.50238 


3.39552 


4.13728 


7.67793 


6.86991 


7.10602 


4.19919 


4.16677 


4.34534 


3(z) 


-6.81708 


-7.11171 


6.87579 


-7.15181 


-4.25387 


4.12951 


4.35574 


-7.33086 


-7.16307 


-8.39439 


4(2:) 


4.97057 


4.98998 


6.55082 


6.15942 


-0.46391 


-0.85430 


-0.68414 


-1.76564 


5.29138 


5.41968 


4(2/) 


4.97057 


4.98998 


6.55082 


6.15942 


2.20720 


2.05072 


2.06772 


1.76564 


5.29138 


5.41968 


4(z) 


4.97057 


4.98998 


6.55082 


6.15942 


9.84219 


9.58102 


9.29718 


11.08793 


-8.87922 


-9.17398 


5(a:) 




10.71204 


8.16900 




-5.11695 


5.14767 


5.30305 


6.25738 


6.31190 


6.21939 


5(y) 




0.00000 


0.00000 




-5.11695 


5.14767 


5.30305 


6.25738 


6.31190 


6.21939 


5(z) 




0.00000 


0.00000 




-5.11695 


5.14767 


5.30305 


6.25738 


6.31190 


6.21939 


6(x) 












-5.83946 


-5.61324 




-8.47946 


-8.50096 


6(y) 












-5.83946 


-5.61324 




-8.47946 


-8.50096 


6(«) 












-5.83946 


-5.61324 




-8.47946 


-8.50096 
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TABLE II. Structural, energetic and electronic features of SRIG 
models: Cross section Q, = 7r6^j„, ratio Rc of concave triangles to 
total number of surface triangles, cohesive energy per atom Ec, 
and HOMO-LUMO gap Eg, [for the clusters that have no gap, the 
degeneracy and filling of the partially occupied level at the Fermi 
energy arc given in brackets]. 



Cluster 




Rc 


Ec (eV/atom) 


Eg (eV) 


33 


52.2 


/ 960 ( %) 


3.816 


0.35 


39A 


100.9 


36 / 1034 (3.5 %) 


3.741 


0.64 


39B 


87.0 


48 / 1156 (4.2 %) 


3.599 


0.77 


45A 


77.6 


/ 1544 ( %) 


3.844 


0.30 


45B 


90.3 


176 / 2896 (6.1 %) 


3.807 


[2 (1/2)] 


49A 


86.0 


36 / 1868 (1.9 %) 


3.858 


0.07 


49B 


81.0 


36 / 1868 (1.9 %) 


3.855 


[3 (2/3)] 


57 


114.0 


288 / 2672 (10.8 %) 


3.712 


0.04 


61A 


141.4 


240 / 2676 (9.0 %) 


3.800 


[3 (2/3)] 


61B 


142.1 


72 / 2128 (3.4 %) 


3.674 


[3 (1/3)] 
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FIGURE CAPTIONS 



FIG. 1: Projected structure of the exterior atoms in the 33 model. Three types of atoms arc 
shown, white apex atoms, hatched atoms and grey atoms. The atoms of each kind are equivalent 
by symmetry. The four triangles outlined by dotted lines represent the four sides of a tetrahedron, 
unfolded around the central triangle. The darker solid lines represent covalent bonds between atoms. 
The five tetrahedrally bonded interior atoms are not shown. 

FIG. 2: Same as in previous figure, for the model 39A. There are four types or atoms, white, 
hatched, grey and black. 

FIG. 3: Same as in previous figure, for the model 39B. 

FIG. 4: Same as in previous figure, for the model 45A. Three types of atoms are shown, white, 
hatched and grey. 

FIG. 5: Same as in previous figure, for the model 45B. The chirality of this model is evident in 
the zig-zag chains of atoms that straddle the dotted lines. This model also serves as the shell for 
the 49 A and 49B models, in which the four apex white atoms or the four face-center white atoms 
become equivalent to adatoms on the (111) (7 x 7) reconstruction, by the addition of a five interior 
atom core. 

FIG. 6: Same as in previous figure, for the 57 model. 

FIG. 7: Same as in previous figure, for the 61A model. Five types of atoms are shown, white 
apex, white face-center, hatched, grey and black. 

FIG. 8: Same as in previous figure, for the 61B model. 

FIG. 9: Density of states of the different SRIG models (shifted on the vertical axis for clarity). 
The height of lines represents the degeneracy (1,2 or 3). The vertical dashed line is the Fermi energy. 
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